Forced Tunneling and Turning State Explosion in Pure Yang-Mills Theory 
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We consider forced tunneling in QCD, described semiclassically by instanton-antiinstanton field 
configurations. By separating topologically different minima we obtain details of the effective po- 
tential and study the turning states, which are similar to the sphaleron solution in electroweak 
theory. These states are alternatively derived as minima of the energy under the constraints of 
fixed size and Chern-Simons number. We study, both analytically and numerically, the subsequent 
evolution of such states by solving the classical Yang-Mills equations in real time, and find that the 
gauge field strength is quickly localized into an expanding shell of radiating gluons. The relevance 
to high-energy collisions of hadrons and nuclei is briefly discussed. 



I. INTRODUCTION 
A. Instanton-Induced Scattering in QCD 

The existence of topologically distinct non-abelian 
gauge fields, with tunneling between corresponding clas- 
sical vacua described semiclassically by instantons j| , is 
one of the most spectacular nonperturbative effects of 
field theory. Significant progress has been made in under- 
standing instanton-induced effects in Quantum Chromo- 
dynamics (QCD), explaining both explicit Ua(1) chiral 
symmetry breaking at the single-instanton level and 
spontaneous SU(iV/) chiral symmetry breaking by the 
instanton ensemble H. Euclidean correlation functions, 
studied phenomenologically and on the lattice, have been 
explained to a significant extent by instantons as well [Q . 

With tunneling phenomena apparently so important in 
virtual quark and gluon propagation, it is reasonable to 
think them also relevant in real processes such as scatter- 
ing or particle production in Minkowski space. We thus 
seek contributions to parton scattering amplitudes from 
the theory of instanton-related objects, and supporting 
experimental evidence. 

With this as our motivation, we concentrate in this 
paper on the theoretical basis of such effects from pure 
Yang-Mills theory. Specific applications to high-energy 
processes with hadrons or nuclei are left for papers to 
follow, although we will discuss phenomenological gener- 
alities where relevant. 

Progress in understanding of the role of tunneling in 
high energy processes has been tempered by technical 
problems for years. Significant insights were obtained in 
the 1980's H and further developed in the early 1990's 
through work in electroweak theory. In this case, 
the instanton-induced cross section is readily identified 
by baryon number violation and many noteworthy fea- 
tures of these processes were found. However, quantita- 
tive estimates of the associated cross sections proved to 
be far below observable limits and interest quickly waned. 
Similar ideas have also been developed in QCD ||, no- 
tably the search for hard processes induced by small-sized 



instantons which continues at HERA |9| . 

Another role for instanton-induced processes has re- 
cently been proposed by Kharzeev, Kovchegov, and Levin 
]l0| and Nowak, Shuryak, and Zahed jy]. These works 
focus on typical QCD instantons, of size p ~ 1/3 fm 
||, which determine the semi-hard scale of Q ~ 1 — 2 
GeV. It was proposed that topological tunneling is be- 
hind the well-known features of high energy scatter- 
ing described phenomenologically by the so-called "soft" 
pomeron. These ideas were further tested in Ref. [p"2| , 
where they were demonstrated to be reasonably consis- 
tent with experimental data. 

Since the 1960's attempts have been made to explain 
high-energy hadronic collisions with multi-peripheral 
models, with various ladder diagrams describing hadron 
production. It was realized that in order to get cross- 
sections which are not falling at high energies, one needed 
vector field exchange in the ^-channel. With the dis- 
covery of QCD, gluons naturally play this role. Generic 
pQCD-inspired models appeared with processes like that 
shown in Fig. |l](a). Eventually this development led to 
the BFKL gluon ladder [[l3| , which produces an (approx- 
imately) supercritical pomeron, a "hard" pomeron with 
the intercept well above 1. Recent studies of high en- 
ergy hard processes, especially at HERA, have indeed 
found strong growth of the cross section with energy for 
truly hard processes (Q 2 3> 1 GeV 2 ), consistent with the 
BFKL treatment. 

But various data at the semi-hard scale of Q 2 ~ 1 GeV 2 
demonstrate rather different growth with energy, consis- 
tent with a "soft" pomeron. Whatever it might be, the 
pomeron should be an object of a particular size deduced 
from the slope of its Regge trajectory, a' ~ 1/(2 GeV) 2 . 
This size of course cannot be explained by basically scale- 
invariant pQCD, and thus calls for a nonperturbative 
derivation. 

Existing models for the soft pomeron also include lad- 
ders made of i-channel gluons, and the differences be- 
tween them lie mainly in the construction of their rungs. 
Each of the various models has a unique answer for what 
is actually produced in gluon-gluon partonic collisions. 
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For example, in Ref. |15[ a pair of pions in the scalar chan- 
nel or a scalar glueball is produced. The introduction 
into this problem of instanton- induced vertices , 
shown schematically in Fig. |l|(b), led to a different idea: 
the object produced is neither a gluon (as in BFKL) nor 
any colorless hadronic state, but rather a colored cluster 
of the gluon field, which in turn decays into several glu- 
ons. It has been shown that the cross section peaks at an 
invariant cluster mass in the range 2.5 — 3 GeV ]lO| , |TT] | . 
It is very important that the states which are produced 
are not a random group of gluons, but rather their co- 
herent superposition. Understanding their composition 
is the main objective of this work. 

A quantum-mechanical interpretation of the collision 
process is central to this question of prompt gluon pro- 
duction. An impressive body of work has addressed this 
problem with classical Weizsacker- Williams fields of glu- 
ons, the Color Glass Condensate jli). Here we con- 
sider a different classical process, one involving topolog- 
ical objects. In Fig. [l](c) we schematically show a bar- 
rier separating two topologically distinct classical vacua, 
with Chern-Simons numbers^] Ncs — and 1. Unlike a 
standard instanton transition, shown by the horizontal 
dashed line, in a high energy collision a finite amount of 
energy is absorbed. This can be viewed as a "forced tun- 
neling" event (either of the other two dashed lines) which 
ends at a turning state, where the total energy is equal to 
the potential energy, so that the paths can exit the (Eu- 
clidean) domain below the barrier. These colored unsta- 
ble objects are close relatives of electroweak sphalerons 
which are defined at the barrier's peak. We will 
demonstrate how these objects then evolve with con- 
served energy, developing into an exploding shell of color 
field. This part of the process is diagrammed with the 
horizontal lines in Fig. [j](c). 

Before we come to these explosions, we will discuss in 
detail the instanton-anti-instanton (II ) configurations 
which describe this forced tunneling. They provide one 
way toward the understanding of the effective potential 
separating topologically different gauge fields, as well as 
the turning states themselves. We then proceed to an- 
other derivation of the same results as static solutions in 
classical Yang-Mills theory constrained in size. The real- 
time decay of the static configuration is studied in detail, 
using both analytic and numerical methods, ultimately 
leading to a description of the expanding shells in terms 
of gluonic quanta. 
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FIG. 1. The top of the figure compares (a) a typical 
inelastic perturbative process (two t-channel gluons collide, 
producing a pair of gluons) to (b) a nonperturbative inelastic 
process, incorporating collisions of a few t-channel gluons with 
the instanton (the shaded circle), resulting in multi-gluon pro- 
duction. The bottom figure (c) shows the same process, but 
in a quantum mechanical way. The energy of Yang-Mills field 
versus the Chern-Simons number, N ca , is a periodic function, 
with zeros at integer points. The instanton (shown by the low- 
est dashed line) is a transition between such points. However 
if some nonzero energy is deposited into the process during 
transition, the virtual path (the dashed line) leads to a turn- 
ing states, from which starts the real time motion outside the 
barrier (shown by horizontal solid lines). The maximal cross 
section corresponds to the transition to the top of the barrier, 
called a sphaleron. 

B. Spherically Symmetric Yang-Mills Fields 

For the SU(2) color subgroup in which we are inter- 
ested, spherically symmetric field configurations of the 
gauge field can be expressed through the following 
four space-time (0, j = 1..3) and color (a = 1..3) struc- 
tures 



*This will be introduced formally below. Here it is suffi- ~ M r ^)®j + B(r,t)Uj + C(r,t)S° 

cient to note only that we consider a definite pair of gauge 
potentials, separated on one of the many coordinates of our 
quantum system. 



A a =D(r,t)- (1) 



with 



e« = ^-, n« = ^-^i, e« = ^. (2) 
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It is convenient to express the scalar functions in Eq. ([j]) 
in terms of four r and t dependent functions, which 
are similar to the fields of the 1+1 dimensional Abelian 
gauge-Higgs model (^=0,1, 4>, a) on a hyperboloid |jl9| : 



A = 



)sma 



(3) 



One can express the field strengths in these terms as 
q = QSj = -[<9 o 0sina + 0cosa(<9 o a - A )]Q? 

+ -[d cj)cosa — (j>sma(d a — A )]Hj 
+ (d A 1 - diA )5L° (4) 



and 



+ -[9i0sina + <j)cosa(dia — j4i)]II° 



1-0 2 



(5) 



where <9 = d t and d\ = d r . 

The action in 3+1 dimensional Minkowski 
space reduces as 



S = / d 3 xdt 
ig 2 



ft)* -ft)* 



47r / drdt 



(1-0 2 ) 2 r 2 2 
— (d A, -dtAo) 



(6) 



with the summation now over the 1+1 dimensional 
(— , +) metric. 

The spherical ansatz is preserved by a set of gauge 
transformations generated by unitary matrices of the 
type 



1 2r 



(7) 



These transformations naturally coincide with the gauge 
symmetry of the corresponding abelian Higgs model: 



a' = a + 0, A' = A p + . 



(8) 



This freedom can be used to gauge out, for example, the 
Ao component. 

Topological properties of the gauge field are governed 
by the topological current 



K — — 
11 ~ 32tt 2 



! ^ vp ° ( Q a vp A a a - 9 -e ahc A a v A h p A c a ) . ( «) ) 



Although this current is not gauge invariant, its change is 
related to the (gauge invariant) local topological charge 



32?r 



2 ^liv^^iv ■ 



(10) 



Within the spherical ansatz and the ^4o = gauge the 
topological current takes a simpler form, 

K° = [(1 - 4> 2 )(dia - A x ) -d x {a- 0cosa)] 



J\ : = [(1 - (f> 2 )d a - d (a - </>cosa)] , (11) 

iharge becomes 
{-daKl-^fra-Ax)] 



87r 2 r 3 

while the topological charge becomes 
1 



+ di [{l~<f){d a a~A Q )\} 



(12) 



Note that only gauge-invariant combinations of field 
derivatives appear here. 

As a "topological coordinate" marking the tunneling 
paths and the turning states one can use the Chern- 
Simons number 



N, 



cs 



d 3 xK = -^ I dr(l-<^H + n - .1,1 



1 x { \\T — OO 

2tt ( a_cosa ^=o 



(13) 



The first, gauge-invariant term is sometimes called the 
corrected or true Chern-Simons number pp| , |2l[ |, Ncsj 
while the second (gauge-dependent) term is referred to 
as the winding number. It is the change in Ncs which 
is equivalent to the integral over the local topological 
charge. 



II. INSTANTON-ANTIINSTANTON 
CONFIGURATIONS 

A. Forced Tunneling 

A brief introduction to the quantum mechanics of glu- 
ons in high energy collisions has been given in the intro- 
duction. The effect of colliding partons can be included 
in various forms. For example, these fields can be rep- 
resented as non-zero external currents which affect the 
tunneling paths of Yang-Mills field. In the zero-current, 
vacuum case, the usual instanton solutions are spheri- 
cally symmetric in four Euclidean dimensions. The colli- 
sion problem of two (or more) partons, on the contrary, 
at non-zero impact parameters does not have even an 
axial symmetry. The reader therefore may wonder why 
this (and all previous works) on the subject consider 3+1 
dimensional spherically symmetric fields. 

The justification for this ansatz is that the absolute 
magnitude of the tunneling field is large compared to ex- 
ternal forces. Also, as will be shown below, spherically 
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symmetric clusters are an energy minimum for fixed size 
and topological coordinate. Should the resulting cluster 
not have exact spherical symmetry one can always ap- 
proach the problem perturbatively, considering first the 
external forces projected onto the direction of tunneling, 
and then other components as small corrections. The 
resulting 1+1 dimensional problem is readily solved nu- 
merically and, to a great extent, analytically. 

Unlike separated instantons (I) and antiinstantons (J), 
combined II configurations are neither selfdual nor anti- 
sclfdual and do not satisfy classical equations of motion. 
They are not extrema of the action, since they describe 
the valley stretching between true extrema - the zero 
field (equivalent to an // at zero separation) and well- 
separated //pair. Substituting any //trial function into 
the of Yang-Mills equation of motion, we find a finite 



(14) 



This means some external current must be applied to 
the gauge fields if we want to use semiclassical analysis. 
The process can only then be interpreted as a classical 
// , or a forced path. There are two interpretations of 
//configurations with different consequences. 

The historical view is that such fields describe quantum 
fluctuations in the Yang-Mills vacuum, the process in 
which a virtual path goes under the barrier, then reverses 
course and ends up in the same minimum from which it 
started. This process has zero net topological charge. 
Naturally, the early studies concentrated on the action 
corresponding to these configurations, the quantity which 
controls its weight in the path integral. The first such 
work was done long ago by Callan, Dashen, and Gross 
|p2| , resulting in a dipole force and the action SS ~ 1/T 4 
at large distance T between the centers. Higher terms 
in the multipole expansion have been discussed in litera- 
ture after that, e.g. [Q. When it was eventually realized 
that quark-induced //pairings are more important for 
the instanton ensemble in QCD [Q, interest in the pure 
Yang-Mills theory waned. 

In this paper we will however take a different view 
of //configurations. Since the external forces from the 
partonic current do work on the //pair, the energy at 
intermediate times is non-zero. We will consider only 
cases in which the fields at positive and negative times 
are essentially the same (modulo a sign and, sometimes, 
a gauge transformation). Thus this energy will be even 
under t — > —t, with a natural maximum at t = 0. As 
a result, all quantities which are odd under this trans- 
formation (like the electric field) naturally vanish at this 
instant. The remaining, purely magnetic configuration is 
what we define as the turning state of this path. 

The resulting action corresponds to an excitation prob- 
ability of this turning state created by the external cur- 
rent J, 



P ~ KOlJlturningstate)! 2 



(15) 



Through this mechanism the excitation of // pairs leads 



to the production of real particles, as advertized in the 
Introduction and to be analyzed in the next sections. 



B. Simple //Trial Functions 

We now consider the simplest example of a possible 
turning state, a straightforward //sum ansatz. With 
it, we can demonstrate some basic features, although we 
will find them insufficient for our purposes and move to 
a more complicated ansatz in the next subsection. 

Written the singular gauge, the sum ansatz is: 



+ 



yi(yi + p 2 ) vlivi + p 2 )' 



(16) 



where we assume that both the instanton and the an- 
tiinstanton (the first and second terms, respectively) 
have the same color orientation and size p. The vectors 
yi = x — zj and yi = x — zj are the distances from the 
observation point x to the instanton and antiinstanton 
centers. In what follows we assume zj = (T/2, 0,0,0) 
and Zj = (— T/2, 0, 0, 0), where the imaginary time be- 
tween centers is T . 

Note that although a single instanton's profile behaves 
as 1/x near the origin, the physical quantity fe) is 
finite. However, for the sum ansatz this feature is lost 
and the same quantity goes as 1 /x 2 near the origin. 

This unphysical feature can be quickly remedied by the 
ratio ansatz fl25| , which for identical sizes and orientations 



S_ \rat: 



>(x) = 



<\a,,jivy\ ^7 



"2 



1 



(17) 



These trial functions are simple enough to have analytic 
expressions for the field strength, the energy of static 
turning states, and the Chern-Simons number. For refer- 
ence, one has the following expressions for the magnetic 
and electric fields squared: 

B 2 = 16384(768i 8 + 1024r 2 * 6 + 3072i 6 + 2304t 6 /? 2 



+6400r 2 i 4 /? 2 



2048r 2 i 4 



h\2r 4 t 4 



1824< 4 /? 4 



4608t 4 /? 2 



l92t 2 R 4 



1024r 6 i 2 

1024r 4 i 2 + lUt 2 R 6 + \2\%r 2 t 2 R 4 



+3072t 4 
+512rW 

+28l6r 4 t 2 R 2 + 288r 4 R 4 + 768r 8 
+A8r 2 R 6 + 3/? 8 )/(16r 4 + 32r 2 t 2 + 8r 2 R 2 + 16t 4 
-8t 2 R 2 + R 4 + 32r 2 +32t 2 + 8R 2 ) 4 , (18) 



768r 6 /? 2 



E 2 = 1048576^ (32r^ 4 + A&t^R 1 + 64r 4 ^ + 64r'H 

12R A + 32r 2 



„2,2 



-80r 2 t 2 R 2 



A8t 2 R 2 



24i 2 /? 4 



-32r b 



Mr 2 R 2 



26r 2 R 4 



-Mr 4 + 12R 2 + 3/? 6 )/(16r 4 



64r 4 /? 2 

+ 32r 2 t 2 + 8r 2 R 2 



-16t 4 - 8t 2 R 2 



R 4 + 32r z + 32t z + 8R Z 



(19) 
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Their scalar product is 



B-e 



-393216ii?(i? 2 + 2 
-2At 2 R 2 + 32r 2 t 2 - 
-8r 2 R 2 )/(16r 4 



4r z 
+ 32t 2 
32r 2 t 2 - 



-8t 2 R 2 



R 1 + 32r A + 32r 



4i 2 )(16t 4 
R 4 + 16r 4 

16t 4 
8R 2 ) 4 , 



8r 2 R 2 



(20) 



where we have set p = 1 and R — T is the intercenter 
distance. 

One can see that, in the simplest case of identical sizes 
and orientations for the I and I, time reflection symmetry 
t — ► — i of the problem is indeed manifest, so that 



A a (r,t = Q) = Q, = Q)=0. 



(21) 



This is illustrated in Fig. g(b). Since configurations of 
this type interpolate between a mostly dual region, with 
E^Zi) — B^zj), to an anti-dual region, where £„(•£/) = 
— B^(zj), it is intuitive that the electric field vanishes in 
the center. 





x,y,z 
T 










) — 


\^_^y t 




FIG. 2. Instanton-antiinstanton configurations, (a) A 
schematic picture in Euclidean space-time. The thick vertical 
line, t = 0, corresponds to the location of the turning state. 
The definition of the inter-center distance T is also shown, 
(b) Distribution along the time axis of 2B%2£% and 2B-S for 
the ratio ansatz with T — p, shown by the solid, dashed, and 
short-dashed lines respectively. The curve for B ■ £ is the only 
one which is t-odd. 



This situation can be readily interpreted in the Aq = 
gauge, in which the electric field is simply the time deriva- 
tive of the gauge field - the canonical momentum in 
Yang-Mills field quantization. Thus the t = magnetic 
state is indeed identified as a turning state, in which mo- 
tion is momentarily stopped. For separation T compa- 
rable to the size p the energy is finite, with a maximum 
E~l/(gp). 

The energy E and Chern-Simons number Ncs f° r ei- 
ther the sum or ratio ansatz can be calculated as a func- 
tion of separation T directly, with the hope that a para- 
metric plot of E(Ncs) w iU reveal a useful profile of the 
barrier as a function of this topological coordinate. 

Alas, for the sum ansatz this idea produces reasonable 
results only for very large separation, T > 2p. When T 
is of the order p, the energy E(T) of the turning state (as 
well as the action for the entire configuration) becomes 
very large, while the topological coordinate Ncs(T) re- 
mains fixed. It is therefore obvious that this set of paths 
does not describe the travel across the ridge separating 
classical vacua which we want to study. Instead, this path 
rises with the barrier but continues to increase as the 
origin is approached, following a direction apparently or- 
thogonal to the topological coordinate we want to study. 

The ratio ansatz yields somewhat better results, with 
finite (and even simple) field structure at all T, including 
the point T = 0. However the results, shown in Fig. ||, in- 
dicate that this set of trial functions can only accomplish 
about one third of the journey we would like to make, in 
terms of the topological quantity Ncs- This inadequacy 
will become apparent after comparison with the results 
to follow. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 

FIG. 3. The normalized energy, ER, versus the 

Chern-Simons number for the ratio ansatz. 
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C. The Yung Ansatz, or Going Uphill 

As a natural set of //configurations, one of us ]26| ] 
suggested starting from a well-separated pair and go- 
ing downhill, along the gradient of the action^]. Nat- 
urally, minimization of the action leads to complete 
// annihilation and zero action. 

It was shown by Yung that these configuration can 
generally be obtained from a solution of the streamline 
equation ^3|. He found solutions for large separation 
T p and used them to derive the next order terms 
in the //interaction, to 0(1/T 6 ). A clever conformal 
symmetry was used to reduce the Yang-Mills problem to 
that of a double well potential. The same trick was then 
used in the numerical solution of the streamline equation 
J^,^7|, in which it was observed that the approximate 
ansatz suggested by Yung also happens to be a very ac- 
curate approximation to true solution, not only at large T 
(as Yung intended) but in fact for all finite // separations 
T. As expected, at T = II annihilation occurs and the 
field strength vanishes^. 

Since we take a different view of //configurations in 
this work, we interpret a solution of the streamline equa- 
tion (or Yung ansatz) as a set of forced paths going uphill 
against the gradient of the force. This process reaches its 
turning point (or state), with some maximal energy and 
Chern-Simons number, and then turns back. Because the 
process proceeds uphill, unlike with other trial functions 
with some arbitrary driving force, we expect that all tra- 
jectories rise along the same path, although those with 
larger T go further up. 

The Yung ansatz for the field configuration is rather 
complicated, and is best written in matrix form: 



igAl m9 (x)=igA Y ^ n9 (x)- 



V2 R (^m - yj )p\ R m 
VmVR 2 vlivi +pT) 
. {^iiVi-v%)p\ . P1P2 



vl + pi 



(o-fiVi ~ Vi) 
R 



V2 



1)2 

R 



1)2 



R 2 V2/2 



(22) 



t This can easily be done numerically, and a set of such 
curves for the quantum-mechanical double well potential and 
the corresponding set of //configurations was found in that 
work. 

* This is not obvious from the Yung expression; it was first 
found numerically. The Yung formula's complicated result at 
T = is nothing but a pure gauge. 



where z is related to the conformal-invariant distance, 
(R 2 + p\ + p\)j{pxP2)- In the case pi = p 2 = p, which is 
the only one we need, this relation reads 



R 2 + 2p 2 + y/(R 2 + 2p 2 ) 2 - 4g 

v 



(23) 



All vectors without an indicative index are SU(2) matri- 
ces obtained by their contraction with the vector er M = 
(1,—ir), for example R = x\ — x 2 — R^cr^. An over- 
bar similarly denotes contraction with a — (l,ir). Note 
that barred and unbarred matrices always alternate, in 
all terms; this is because one index of each matrix is dot- 
ted and the other not, in spinor notation. Finally, the 
additional coordinate is 



2/2 = x 2 



RP2 



ZP1 - P2 



(24) 



Note that the first term is the instanton in the sin- 
gular gauge, the second is the anti-instanton in the reg- 
ular gauge, and the third is a "correction" term. The 
benefit of this representation is that the same 't Hooft 
symbol appears in all three terms, and the entire con- 
struction originates from conformal transformation of a 
spherically symmetric configuration in which /, / share 
the same center. An unfortunate feature of this expres- 
sion is that time-reversal symmetry is far from obvious, 
and it is not clear that the electric field at the mid-plane 
vanishes. However, this is in fact the case and the field 
at t — can be interpreted as a turning state. 

Although all three trial functions are similar at large 
//separation T, they are drastically different at T ~ p. 
The Yung ansatz is the only one which allows us to rea- 
sonably study the effects of a large change in topological 
number. The variation of the Chern-Simons number of 
the turning state [t = 0) as a function of the // separation 
T can be seen in Fig. H. In this case we scan the entire 
range [0,1]. 




FIG. 4. Ncs versus the distance between II centers T in 
the Yung ansatz. 
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We now proceed with a more detailed study of the 
static turning states, residing on the t = 3-plane. The 
simplest observable is the shape of the corresponding 
magnetic field squared, or the energy density distribu- 
tion, shown in Fig. || for few selected values of II dis- 
tance T. Note that the curve for T = 2 (the most like 
the sphaleron) show indeed the largest magnitude of the 
magnetic held. The shape is however rather uniform. 
Note also that, unlike the case of the faulty sum and ratio 
trial functions, for smaller T the field strength decreases, 
ultimately disappearing at T = 0. 



Energy*Size versus Ncs 




FIG. 5. The B(r) 2 profile, not normalized, for the four 
values of the II distance T (in units of p) indicated in the 
legend. 

The energy and energy density of the turning state 
configurations is therefore rather different for different T. 
However, as seen from Fig. ||, the physical sizes of these 
objects are different as well. As classic Yang-Mills theory 
has scale invariance, one may wish to make the more 
natural comparison of a scale-invariant combination, the 
energy times the r.m.s. radius, R, defined as 



R 2 



Jd 3 rr 2 B 2 



Jd 3 rB 2 ' 

In these terms, the normalized energy is 



ER= - 
2 



d 3 rr 2 B 2 x / d 3 rB 2 



1/2 



(25) 



(26) 



This quantity is plotted versus the topological charge dif- 
ference in Fig. ^, and indeed displays a parabolic-looking 
maximum near Ncs = 1/2. 




Time dependent paths: E*R vs Ncs 




FIG. 6. The normalized energy, ER, versus the 

Chern-Simons number for the Yung ansatz. Plot (a) shows 
the positions of the turning states for various T, while (b) 
combines many points along the path (t 7^ 0); their small 
spread means that Yung ansatz is nearly going directly up- 
hill, thus passing via the same points for different T. 

Instead of only looking at the static t = (and zero 
electric field) turning states, one can instead follow the 
(scale invariant) energy ER and the Chern-Simons num- 
ber as a function of time t along each each path. As 
expected, all the paths in Fig. ^(b), for any T, actually 
climb nearly exactly the same cliff, as they propagate into 
larger values of our topological coordinate. 



III. TURNING STATES FROM CONSTRAINED 
MINIMIZATION 

We will now define turning states in terms of the gauge 
field, which connect the Euclidean and Minkowski do- 
mains of the field's path. The turning state is character- 
ized by the condition that the generalized momentum, 
which in the Ao — gauge coincides with the chromo- 
electric field, vanishes or, equivalently, that all first time 
derivatives of the spatial field components are zero. Us- 
ing the notation introduced in Section IB, this in turn 
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means that do4> = d$a = OqAi = at the time when the 
transition occurs. From now on we assume that moment 
to be t = 0. 

At any given time it is possible to use the special gauge 
transformation, Eq. um, with a time- independent angle 
j3 to gauge out Ai(r), still within the A — gauge. At 
t = the energy of the field can thus be written 



E 



47T 



dr 



{d r cj)) 2 +(j) 2 (d r a) 



(W 2 ) 2 

2r 2 



(27) 



We now address the question of the minimal potential 
energy of static Yang-Mills field, consistent with the ap- 
propriate constraints: (i) a fixed value of the (corrected) 
Chern-Simons number, and (ii) a given value of the r.m.s. 
size. The former parametrizes the position of configura- 
tion on the topological scale, and the latter is needed to 
break dilatation symmetry of the problem, which other- 
wise prevents any configuration of finite size from being 
the minimum of the energy. 

We will break the scale invariance of the theory by 
setting a requirement that the ratio 



J d 3 xr 2 B 2 
J d 3 x B 2 



(28) 



has a particular value, p 2 , for the static solution we seek. 
To keep both the Chern-Simons number and mean radius 
constant we introduce Lagrange multipliers and search 
for their minimal combination ofpl 



The Euler-Lagrange equations for the remaining fields 
are 



dU - <X^a) 2 



(l-0) 2 
cos 2 £ 



+ 2K<j>d t a = (31) 



and 



d^ 2 d ia ) + Kd^l - (f> 2 ) = . 



(32) 



Finiteness of the energy requires the boundary conditions 
(f> 2 (-n/2) = </> 2 (tt/2) = 1. Eq. © integrates to 



df a = — k 



1 



(33) 



with a vanishing integration constant as follows from the 
form of the energy. After the substitution of d^a into 
Eq. p|) one has 



5 cos z £ 



.1 



(34) 



A solution to this equation exists for — 1 < k < 1, 

4> 2 = 1 - (1 - k 2 ) cos 2 £ . (35) 

Hereafter we assume that (f> is positive. 

In term of the usual r coordinate, we have instead 



l-(l-K 



2 . 4p 2 r 2 N 1/2 



(r 2 + p 2 y 



dr 



r 

1 + - 

P 2 



(d^) 2 + cj) 2 {d r a) 2 + 



(1-0 2 ) 2 

2r 2 



>1 



/ dr(l ~ (f> 2 )d r a. 



(29) 



where the tilde denotes the constrained energy. It is con- 
venient to introduce a dimensionless variable, 



£ = 2 arctan 



2 " 



so that 



E : 



8tt 



tt/2 



di 



-tt/2 



+k(1 - (f> 2 )d^a 



where k = rjpg 2 /(32ir 2 ). 



(30) 



d r a(r) — —2k- 



4> 2 r 2 + p 2 

The sphaleron solution corresponds to n = and 



(36) 



0(0= K, P l l ■ <*{r)=-K6{r-p). (37) 

For any k and mean squared radius (r 2 ) = p 2 , the 
potential energy density is 



~ - 24P {r 2 +p 2 Y 



(38) 



the integral of which is the potential (magnetic) energy 
of the static configuration, 



E B = 3tt 



2 (I -* 2 ) 2 
9 2 P 



(39) 



The corrected Chern-Simons number, computed from the 
first term of Eq. (0) , is 



" Without the term introduced to fix the Chern-Simons num- 
ber or, equivalently, for zero corresponding Lagrange multi- 
plier, this problem would be equivalent to the SU(2) sphaleron 
on a 3-d sphere that was solved by Smilga feg|. 



A> C s = i sign( K )(2 + | K |)(l-| K |) 2 . 



(40) 



Figure ^ shows the profile of the potential energy Eb 
versus Ncs- It is very similar, although not identical, to 
the findings of the preceding section (see Fig. ^) where 
Yung's ansatz was used for forced paths. 



8 




FIG. 7. The potential energy Eb versus Ncs, for the 
analytic turning state solution of Eq. (jirj ) . 



IV. EXPLOSIONS OF THE TURNING STATES: 
ANALYTIC TREATMENT 

We are now going to use the static field configuration, 
found in previous section, as an initial condition for real- 
time, Minkowski evolution of the gauge field. Let us first 
consider the equations of motion in the 1+1 dimensional 
dynamical system. Variation of the action, Eq. (0) , gives 



(1- 







[d^a-AJ] =0 



<p 2 (d ia -A 1 )-d 



<j ) 2 {doa-A )-d 1 



-(d A 1 -d 1 A ) 



-(d A 1 -d 1 A ) 



= 
= 0. 



(41) 
(42) 

(43) 



The solution of Eq. (42) has the form 



4> 2 (d a-A a ) 
c)j 2 (d 1 a-A 1 ) 



(44) 



where i/j(r,t) is an arbitrary smooth function. Eqs. ( [43] ) 
are consistent with this solution if 



d A 1 -d 1 A = -^ 



Now, combining Eq. 



and Eqs. (O) one has 



d A 1 - diAo 



2i/j 



(45) 



(46) 



which can be viewed as a necessary and sufficient con- 
dition for f/j to be a solution for Eq. p^ ) and Eqs. ( f43| ) 
simultaneously. Eq. ( |4l| ) is now 



(d^) 2 , (W 2 



= 0. 



(47) 



The initial conditions for Eqs 
(r), 



j) and (47) are 

4>{r,0) 
d (t)(r,t)\ t=0 = 0, 

d r tp(r, 0) = -4>{rfd a{r) = => ip(r, 0) = 0, 
5o^(r,i)| t=o = -0(r) 2 9ia(r), 

where the ^independent fields on the right sides of the 
equations are the static solutions of <f> and a from the 
previous section. 

As with static solutions, it is more convenient to dis- 
cuss the time-evolution equations in hyperbolic coordi- 
nates. Let us choose w and r such that 



p cos lo 



cos r — sin lo 



t 



psmr 



(48) 



The physical domain of < r < oo and — oo < t < oo is 
covered by — ir/2 < lo < ir/2 and — it/2+uj < t < tt/2—lo. 
For t > 0, the corresponding domain is — ir/2 < lo < tt/2 
and < t < tt/2 — lo. This change of variables ( [48| ) is a 
conformal one. 

In the new variables Eqs. (^) and ( p7| ) become^] 



K4> + K4> 



a 



2^ 







= 0. (49) 



Before solving these equations let us note that it is 
possible to predict the large-i behavior of gauge field from 
the form of the conformal transformation ([18]) . Indeed, 
the t — > oo limit corresponds to the line t — n/2 — lo on 
the (lo,t) plane. If one now takes the limit \r — t\ — > oo 
(regardless of the limit for \r — t\/t), the position on (lo, t) 
plane is either lo — ► — n/2 , r —> or lo — > tt/2 , r — > n. 
This means that the entire line r = 7r/2 — lo corresponds 
to space-time points with finite differences between r and 
t and, therefore, if <j> and ip are smooth functions of lo and 
r, then for asymptotic times the field is concentrated near 
the r = t line. This corresponds to the fields expanding 
as a thin shell in space. 

We must now supply Eqs. ( (49| ) with initial conditions, 
which are 



<P(lo, t = 0) 2 = 1 - (1 - k 2 ) cos 2 LO 
9 T 0(w,r)| r=o = 
i/)(u,t = 0) = 

d T ip(Lo,T)\ r=0 = - — d t <tp(uj,T)\t=o 

1 — sin lo 

= k(1 — k 2 ) cos 2 LO . 



(50) 



** One can find a discussion of Eqs. 
solutions in IE if. 



and some of its 
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One of the solutions of Eqs. (fpfh first found in 1977 by 
Liischer [E9|| and Schechter [(301, is 



1 - (1 - <7 2 (t))cos 2 uj 



cos 2 UJ . 



with a function q(r) that satisfies 

q-2q(\-q 2 ) =0. 



(51) 



(52) 



This is the equation for a one-dimensional particle mov- 
ing in double-well potential of the form U(q) — (1 — 
<Z 2 ) 2 /2. 

We now have to check that the Luscher-Schechter so- 
lution satisfies the initial conditions, (|50|). This is indeed 
the case if one identifies q(0) — k and takes q(0) = 0. 
For the initial condition of this type (i.e. for energy 
e = q 2 /2 + U(q) < 1/2), the solution of Eq. (||) is 



q(r) = qdn(q(r - T ),k) , 



(53) 



where dn is Jacobi's function and q — y/2 — k 2 is the 
second stopping point for a particle in the potential U (q). 
We have also defined 

,2 1 — ft 2 , „ T 
k = 2 ;r and Tnq = — , 

2 — k m 2 ' 

where T, the period of oscillations in the potential U(q), 
is T = 2K(k), with K(k) being the complete elliptic 
integral of the first kind. The idea is, of course, that 
"oscillations" in r begin from the rest point, close to 
r = 0. 

Let us now look at several properties of the solution 
for large times. The solution ( p3| ) is apparently regular in 
the (uj, t) plane, and therefore for large times the field is 
concentrated near r = t. At asymptotic times the energy 
density, e(r, t), is given by 



47re(r, t) 



9 2 P 2 



[l-n 2 f 



^p 2 + (r - t) 2 
The change in topological charge is 



(54) 



OO 

AQ = J (Pxdtd^K* 



— / drdt 
2tt 



— k(3 — n 2 ) — sign(ft) arccos 



/ cn(g7r, k) 



(55) 



2 v \ v dn(q ; 7r, k) 

The evolution of Ncs begins from time t = 0, where 

N CS (0) = Jsign( K )(l - |«|) 2 (2 + |«|) , (56) 
and as t — > oo its limit is Ncs(°°) = ^cs(0) + AQ. 



We now estimate number of gluons produced by the 
described evolution. In 0, ip language the chromoelectric 
and chromomagnetic fields are 

p„ 1 (p. , • d r ipcosa\ 
E" = - | dt<psma ; 1 6" 



i(^cosa + ^W3^ (57) 



+ i^sina-M^)n ? + i^S«. (58) 

Terms proportional to E" are longitudinal and die out as 
t — > oo. The remainder is a purely transverse field. Now 
let us take into account that for t — > oo, d r <p — > dt4>, and 
the same for ip. Therefore 

™ 1 fa j. ■ d r ipcosa\ 
Es — ► - drdsma ; 6? 



1 /_ <9 r i/>sina\ 
— I d r (/)cosa + ) nj, 



™ 1 ( r, , <9 r i/>sinoA ^ n 

1 /„ , . cos a \ 

+ -(^ S1 na- y jn°. 



(59) 



(60) 



The main result becomes apparent when we choose a 
gauge where 



%<p cos a + d r ip sin a = . 



in which 



1 — K 



p 



Bl 



rp \p 2 + (r — <) : 



1-K 2 



.3/2 



e?. 



3/2 



rp \p 2 + (r — t) 2 
We now perform a fourier transform, finding 

Ef(js) = 4np{l - K 2 )K 1 {ujp)e^ 
B a Ak) = ^p{l-K 2 )K 1 {ujp)m, 



(61) 



(62) 



(63) 



where 0° and II" are the color/space projectors in mo- 
mentum space analogous to those in coordinate space 
the frequency uj = \k\, and Ki is a Bessel function. One 
can easily verify that — CjimkiE^/k, as is required 
for a radiation field. 
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V. EXPLOSIONS OF THE TURNING STATES: 
NUMERICAL STUDIES 



+ 



In the previous section we studied the asymptotic be- 
havior of turning states constrained by sphaleron size 
and Chern-Simons number. In this section we consider 
step-by-step evolution of the turning states, a numer- 
ical analysis similar to sphaleron decay in electroweak 
theory . The numerical approach naturally allows 

for mathematical flexibility, and it is used to consider 
the decay of static states which replace the unphysical 
power-law behavior of the fields at large distance with 
a phenomenologically more appropriate exponential tail. 
The classical field configurations are thus fixed in size 
indirectly by a mass parameter, a constraint which is 
subsequentially relaxed as the state quickly decays into 
free-streaming gluons. 

As a result of scale invariance, the QCD instanton is of 
indeterminate size. While it is clear from phenomenology 
and lattice studies that the instanton vacuum favors a 
somewhat narrow size distribution centered at p ~ 0.3 
fm, the reason for this is yet unknown, although it is 
presumably due to interactions between instantons. It is 
thus natural that related classical objects, born in some 
way from the excitation of instantons, share a similar size. 
We arrange this by introducing a phenomenological gluon 
mass term in the initial configuration which is promptly 
relaxed as this unstable configuration begins to decay. 
We stress that the relaxation of this size constraint does 
not initiate the explosion; we will show that the turning 
state is an unstable configuration regardless of the mass 
term's presence. 

The similarity between this procedure and electroweak 
sphaleron decay is clear, but not mathematically con- 
tinuous. For the Higgs mechanism, the sphaleron size 
is constrained by the scalar vacuum condensate which 
introduces an effective mass for the gauge fields. This 
condensate vanishes at the sphaleron center, where the 
classical gauge field action is maximal. This feature per- 
sists in the limit of infinite gauge-Higgs coupling, when 
the Higgs field is fixed at its VEV constant for all other 
points in space. Inserting a mass term for the gauge fields 
by hand, as we do here for QCD, is thus very similar to 
this infinite coupling limit of electroweak dynamics, the 
only difference being that the mass is finite in all of space. 
This leads to a difference in field behavior at the origin. 

We begin with the Yang-Mills action, with a phe- 
nomenological mass term added, and will look for static 
solutions with spherical symmetry in Minkowski space. 
The action is written: 



5' 




W) 2 -(^) 2 



l</> 2 ) 



As before, we have taken a spherical ansatz similar to 
Witten's |1!| for the gauge field. Following the notations 
ofEq. (§), 



.4 



D 



, C = A U D = A , (65) 



where all scalar functions depend on both r and t. We 
initially work in the temporal gauge, where Ao(r, t) = 0. 

The equations of motion are easily obtained, and we 
find 

k - <P'I + ^(<i>l + <t>l-l)+ <t>iA\ - 2$Ai 
-<Mi + m 2 (0! + 1) = , 

& " + ^ fal + ^ ~ *) + + 2( t>'l A l 

+^iA[ + m 2 4> 2 = , 

^'1 + (^1 + + ^2 (^i^2 - <h<&) 

+m 2 A 1 =0. (66) 

In the static limit, we seek a purely magnetic turning 
state solution. One can be found for the 4>i field compo- 
nent from the equation 



1 



61 + 1) = 0, 



(67) 



which is simply Eqs. ( p6[ ) with 02 = A\ = and written 
in terms of the dimensionless variable x = mr. With the 
boundary conditions 

(f>i(x) —> 1 as x — > , 
</>i (x) — > — 1 as x — > 00 , 

A numerical solution is easily obtained and shown in 
Fig. ||. This is quite similar to the approximate elec- 
troweak solution found by Klinkhammer and Manton fT§| ] 
in the limit of infinite Higgs self-coupling. The primary 
difference is the behavior near the origin, which in the 
QCD case involves a logarithm for x <C 1: 



4>i(x) — 1 H — x 2 \tlx - 
3 



0{x A \nx) 



(68) 



where a = 1.98 was determined numerically. Defining 
the turning state's size as the radius at which the profile 
is at half its maximum (i.e. where it crosses the origin), 
we find mp ~ 0.9. We match this to the size of the 
average instanton and find 



(64) 



m 



0.9p" 



540 MeV. 



(69) 
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FIG. 8. The static turning state solution, (j>i(r, 0). 

Before we discuss the decay of this configuration we 
must find the unstable modes orthogonal to it which de- 
termine the "downhill" directions in field space. This is 
done by solving eigenvalue equations for fluctuations in 
the fields 4> 2 and Ax in the presence of the turning state 
configuration. 

We take the terms linear in <j> 2 and A\ from Eqs. ( |66| ) 
and require 

fo{x,t) = -u; 2 (t) 2 (x,t), A 1 {x,t) = -uj 2 A 1 (x,t). (70) 
We then have the eigenvalue equations: 
1 ' 



J 2 ' 



1 



I Mi 



b 2 + 2cf)' 1 A 1 +4> 1 A' 1 = 0, 
- -l(^ 2 -^4) = 0) (71) 



where <f>\ is the classical solution in Fig. g and the di- 
mensionless frequency is f2 = uj/m. The longitudinal 
field component may be eliminated with 



A x 



2 - fr^) 

(n 2 ~i)x 2 - 2(t>\ 



(72) 



Substituting this into the first of Eqs. (|7l|), we find the 
behavior near the origin: 



4> 2 {x) =cxi( 1 - VT ^) 



(73) 



where 



V 



i + n 2 
l-n 2 



and c is an arbitrary normalization. Both fields vanish 
as x — > oo. 

We have solved these equations numerically finding 
the wave functions plotted in Fig. ^| with the frequency 



u> 



-3Am 2 



(74) 



The function A\ is logarithmically divergent at the ori- 
gin, reflecting the difference between this massive model 
and the electroweak case Q in which the Higgs always 
vanishes at the origin^. There is therefore no smooth 
continuation between this model and the electroweak in 
the limit of large coupling. 
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FIG. 9. The unstable eigenmodes, 4>i(r) and Ai(r), arbi- 
trarily normalized. 

These solutions for the unstable modes, which along 
with the classical 4>\{r) complete our initial conditions, 
were put on a lattice with spacing Ax = 0.01 and evolved 
at time steps of At = 5 x 10 -4 , where r = mt. The un- 
stable modes, acting as a small push to properly initiate 
the decay, were normalized as 



dx ((/)2(x) 2 + 2x 2 A 1 (x) 2 ) = 5 x 10" 



(75) 



Coincident with this push we set m = 0, in effect turn- 
ing off the mass term, since here we are interested in the 
turning state decaying into the vacuum where no such 
term is motivated. Although this effectively removes the 
size constraint on <f>\{r 1 0), the subsequent dynamical ex- 
pansion is a result of the push in the unstable directions 
of this saddle-point solution rather than an inflation of 
the classical field configuration. 

Once the real-time solutions to Eqs. (66) with m — 



have been found at a given time step, the total energy is 
readily computed as 



E 



-in 
7/2 



dr 



1 2 i 2 
2 



2r 2 



- 2 Ax (<j)i(/>2 - </>i0' 2 ) 



The dominant unstable mode in the electroweak as A — > oo 
is the scalar field component orthogonal to the condensate 
[p3l. This mode is absent here. 
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-A\ 



(76) 



At every time step it was verified that the energy remains 
equal to that of the initial state. Taking the instanton 
vacuum value of 



(77) 



the total energy of the decaying object was calculated 
and found to be 



E = 4.62/r 1 ~ 2.8 GeV. 



(78) 



The Chern-Simons number was also computed at each 
time step. In our present gauge, this is written as 



N, 



cs 



= ^J dr [( 1 - <h) <t>2 + M2 - (1 - <t>\ - 4>l) At] . 

(79) 

The gauge invariance of changes in this quantity were 
verified numerically. 

The energy and Chern-Simons densities, defined as the 
integrands of Eqs. ( f76| ) and (f79|), are shown in Fig. [l(]. 
The shell-like expansion is illustrated in these plots, as 
well as the similarity of the energy profile at all times. 
Note that these plots are in a 1+1 dimensional descrip- 
tion, and the corresponding three-dimensional radial den- 
sity differs by a factor of 1/r 2 . 




-0.05 



FIG. 10. Energy (solid lines) and Chern-Simons number 
densities (dashed lines) for three times during the explosion, 
t =0.4, 3, and 6 fm. (Note the scale difference between the 
two quantities.) 

It is clear from the curves in Fig. [lo] that the Chern- 
Simons number changes during the decay. As shown in 
Fig. [y], our fields stabilize at long times with ANcs — 
0.12. From this we see that the turning state does not 
complete an instanton transition, which would require 
a return to a state with integral Chern-Simons number. 
Due to this freezing in the topology, nontrivial fermionic 
solutions will accompany the resulting Yang-Mills fields. 
These will be discussed elsewhere. 
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FIG. 11. Change in the Chern-Simons number. 

The transition from a purely magnetic configuration 
to one with equal electric and magnetic components is 
shown in Fig. |l2j, for an early and late time in the evolu- 
tion. The decay progresses rather rapidly; at t = 1.4 fm, 
the ratio E(r) 2 /B(r) 2 ~ 0.95 for all r. Thereafter the 
ratio continues to quickly approach unity. 
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FIG. 12. The electric (solid lines) and magnetic (dashed 
lines) fields at times t — 0.4 fm and t = 3 fm. Plotted are 
jX 2 E(x) 2 and ^x 2 B(x) 2 ; their sum is the energy density. 

In order to analyze the final state at late times, we we 
work in a gauge in which B(r, t) = 0. This requires the 
transformation to a new set of fields: 

0i = 4>\ cos 9 + 4>2 sin 9 
4>2 = —<fii sin 9 + cf>2 cos 9 



where 



Ax=A x - 9' 

Ao = -e, 



9 — arctan 



(80) 



(81) 



Promptly dropping the tildes, we write the total energy 
in terms of the new fields, 
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47T , 

E = — I dr 



4>i +K) 2 + 



(i-#) S 



2r 2 



r 

~2 



(i, 



A' 



(82) 



At late times, the field strength is confined to a thin 
shell at radius r = t. Free, expanding field behavior is 
also observed numerically, such that 



$(r,£) = -$'(r,t) 



(83) 



for each field <j>\, A\, and A Q . This simplifies the equa- 
tions of motion for the latter two, 



Ai-Aq = — 2 4>\ Ai 



AX -A! =— A'-A 1 



(84) 



in that the right hand sides of both vanish as the shell 
expands at large times. We can thus conclude that 

ii - A' = (85) 

for mt 3> 1. 

The condition for this gauge is 

Aq4>i + 2A 4>! - A[(j)i - 2A 1 <f>' 1 = . (86) 

Combining this with Eq. (|S5|), we have 

0! (i - 4) - 2 (j)[Ai - faAo) = , (87) 
and can deduce that 

A + Ax = (88) 

at late times. 

The contribution from (f>i will be of the form 

fa(r,t) = l + <p(r,t), (89) 

where ip is an excitation above the vacuum. Using the 
result just obtained (|88|), its linearized equation of motion 
simplifies to 



<P - <p" + ^ = . 



(90) 



The solution is of the form 

tp(r,t) — J dkip(k)r [ji(kr) + iyi(kr)] cos(kt) , (91) 

with fourier amplitudes (f(k) and the spherical Bessel 
functions ji(z) and yi(z). 

Although the majority of gluon radiation is carried in 
the 4>\ field, physical quanta also lie in small excitations 
of the field 



tp(x,r) = 0i(x,t)A o (x,t) , 

which encodes oscillations between the transverse degrees 
of freedom. From Eqs. (|8^) and (|S8|), we have a wave 
equation, 



<$ - = . 



(92) 



These harmonics contribute to the energy via the final 
term in Eq. @. 

The total energy can then be written in momentum 
space as 



E 



16to 



(93) 



in terms of the dimensionless momentum p = k jm. The 
fourier amplitudes are computed from the solutions of 
the spatial fields: 



dxpx [ji{px) + iyi(px)} <p(x) 



dxe ipx tp(x) 



(94) 



Numerically, this expression for the energy is within 1% 
of that of the initial configuration for all times t > 3 fm, 
further demonstrating the rapid onset of free-field behav- 



VI. PRODUCTION OF TURNING STATES IN 
HEAVY ION COLLISIONS 

A. The pQCD Cutoff in Vacuum Versus Excited 
Matter 

Although this paper does not generally deal with phe- 
nomenological applications, we will eventually come to an 
estimate of the number of gluons produced in the explo- 
sion of a turning state. The resulting divergence in this 
number cannot be resolved without some explanation of 
the limited applicability of the classical Yang-Mills de- 
scription. This leads to the issue of the pQCD cutoff. 

It is well known that in the QCD vacuum the "semi- 
hard" or "substructure scale" Q 2 ~ 1 — 2 GeV 2 is simul- 
taneously the lower boundary of pQCD as well as the 
upper boundary of low energy effective approaches such 
as chiral Lagrangians. Furthermore, in any discussion of 
instanton-induced reactions it is implicitly assumed that 
instantons were the primary source of that scale, and 
since they are included explicitly no other nonperturba- 
tive cutoffs are needed. This, of course, is not strictly 
true as confining forces require the final state be com- 
prised of hadrons, but we make the usual separation of 
scales and assume that final state interactions merely re- 
distribute the wave functions without changing the total 
probabilities. 
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We assert, however, that the pQCD cutoff is quite 
different in heavy ion collisions. It has been argued 
over the years that excited matter might be in a Quark- 
Gluon Plasma (QGP) phase of QCD (see e.g. Ref. @). 
Whether equilibrated or not, it is nevertheless qualita- 
tively very different from the QCD vacuum: instantons 
are suppressed, and there is neither confinement nor chi- 
ral symmetry breaking to set a nonperturbative cutoff. 

Therefore the limits on Yang-Mills field description are 
entirely different, and actually determined by much sim- 
pler phenomena. The QGP, a plasma-like phase, screens 
itself perturbatively |$4j . A quasi-particle description be- 
comes appropriate, in which the quarks and gluons have 
finite effective masses. In equilibrium and at high tem- 
perature these are "thermal masses"; the gluon has the 
well-known effective mass 34 



9 2 T 2 



M, 



3 



Nf 
6 



(95) 



where N c and Nf are the number of colors and flavors, 
respectively. Although this mass grows with temperature 
at high T, just above T c it is actually smaller than the 
pQCD cutoff in vacuum. Such non-monotonic behavior 
is confirmed by lattice thermodynamics data, which can 
be well fitted with quasiparticle masses. 

Moreover, in the "RHIC window", T c < T < 3T C , one 
finds the approximately constant gluon and quark effec- 
tive masses p5| 



physical relevance, as the gluonic field is not massless 
phenomenologically (and even vacuum instantons ought 
to have exponential tails). We now consider both results. 

To find the number of gluons in each mode one com- 
pares the field strength in the momentum representation 
to those which have energy to. In the evolution described 
analytically in Section IV, the occupation number is 



v{k) 



64tt 



.9 



(97) 



and the gluon energy distribution function is 
n( W ) = ^(l-K 2 )W^i 2 M 

g 2 

The corresponding energy spectrum E(u) — Lun(uj) is 
shown in Fig. 



M g w .4 GeV , 



M q « .3 GeV, 



(96) 



the first of which provides the cutoff of our classical treat- 
ment; at this scale the classical Yang-Mills action for glu- 
ons is to be modified by inclusion of an appropriate ef- 
fective Lagrangian describing such screening effects, such 
as those suggested by Taylor and Wong j36| . 



B. Multiplicity and Spectra of Prompt Gluons 

With solutions for the fields at all times, in both the 
analytic and numerical treatments of the previous two 
sections, we have analyzed the final states to determine 
the particle number. While we find a similar number of 
produced particles from both approaches, about four glu- 
ons produced per turning state, the energy distribution 
of this prompt glue is very sensitive to the classical con- 
figuration used to initiate the explosion. While the field 
solutions found using constrained quantization (Section 
III) and an effective mass (Section V) are qualitatively 
similar, both leading to an expanding shell of radiation, 
we find a substantial difference in final state momentum 
distributions. This can be traced to the details of the 
two solutions as r — ► oo, where we contrast a power-law 
behavior in Eq. (^6|) with an exponential fall-off in the 
solution of Fig. g, as 4>i(r, 0) ~ exp(— mr) — 1. As men- 
tioned above, we consider the second case to be of greater 



FIG. 13. The energy spectrum E(u>) obtained from the 



analytic solution, in units of pr(l 



■ up. 



Because the fourier transforms of the fields are finite, 
the occupation number (^) behaves as l/to for small 
u). The number of particles is thus logarithmically diver- 
gent and should be cutoff at some low scale where the 
pQCD description of gluons is no longer valid, leading 
to N g ~ log{j^—) . As explained in the previous Subsec- 
tion, in heavy ion collisions this cutoff, identified with the 
gluon effective mass of Eq. (^6|), is still relatively small 
as compared to the typical momenta of gluons produced. 

The prompt particle energy distribution was also ob- 
tained in the numerical treatment of Section V, defined 
as the integrand of the expression in Eq. (|9^), is shown 
in Fig. [[4] at a very late time in the evolution (r = 50 
or t ~ 20 fm). Like the previous distribution, it is finite 
at the origin, but in contrast it peaks at nonzero mo- 
mentum. For illustrative purposes it is compared with a 
thermal distribution of bosons at a temperature T = 285 
MeV. We note that in an equilibrated environment, the 
effective mass and screening effects will only modify this 
profile below k — 0.4 GeV, where a relatively small part 
of the spectrum resides. 
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The produced gluons are free streaming, with no mech- 
anism for equilibration, and yet our distribution is very 
similar to the thermal one for momenta below about 1.5 
GeV. Such a nearly thermal distribution has also been 
obtained from similar calculations using an entirely dif- 
ferent classical field approach in Ref . |57j . One can specu- 
late that such an ostensible equilibration may contribute 
to the the success of hydrodynamics in calculating par- 
ticle spectra and elliptic flow at RHIC |38|]. Our finding 
of four physical gluons from every decayed turning state 
is also in line with RHIC entropy production, assuming 
that the density of these classical objects corresponds to 
that of instantons in the vacuum. The total energy of 
the turning state, found above to be 2.8 GeV, is carried 
by gluons with a distribution peaked around 800 MeV. 
This average can be viewed as an upper bound, since in a 
more complete treatment a portion of this energy will be 
used for the production of light quark pairs. This next 
step will be addressed in a later publication. 



perhaps deconfined medium. Observable consequences in 
both situations were discussed recently in Ref. and 
the results of this work support many of the estimates 
therein. Finally, we have obtained energy spectra for the 
prompt gluons that serve as the initial state in the dy- 
namics of a heavy ion collision and found that, unlike the 
overall explosive dynamics, it is sensitive to the details 
of the initial turning state profile. We plan to further 
investigate phenomenological implications of the turning 
states in later works, with the role of fermion production 
a top priority. 
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FIG. 14. Energy spectrum of prompt gluons (solid line), 
obtained from the numerical solution, and a thermal distri- 
bution with T = 285 MeV (dashed line). 



VII. CONCLUSIONS 

In this work we have studied forced tunneling in 
pure SU(2) Yang-Mills theory. This process, generated 
through the excitation of instantons in the QCD vac- 
uum, leads to unstable classical turning states which ex- 
plosively decay into gluonic radiation. These states and 
their decays are similar to the physics of sphalerons in 
the electroweak theory. 

If this semi-classical treatment of pure Yang-Mills the- 
ory is indeed relevant to the physics of QCD, the turning 
states should play a prominent role in the production of 
glue in high-energy hadronic collisions. In the case of 
NN collisions the produced gluons will propagate into 
the QCD vacuum, to be quickly recombined into sec- 
ondary hadrons. For heavy ion collisions, however, the 
large quantity of prompt glue produced from the many 
turning states would be released into a highly excited, 
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